void drawOBEtaCorr(const Int_t ijb)
{
  const TString jbs[]={"jetPt05_10", "jetPt10_15", "jetPt15_20"};
  ((TDirectoryFile*)(gDirectory->Get(jbs[ijb])))->cd();

  TH3D * hh=(TH3D*)gDirectory->Get("hCorr");

  const Int_t npt = hh->GetNbinsX();
  const Int_t ntype = hh->GetNbinsY();
  const Int_t neta = hh->GetNbinsZ();

  const Int_t itype = 0;//proton

  for(Int_t ieta=4; ieta<=neta; ieta++){


  }
}
